Response maxima in modulated turbulence 
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Isotropic and homogeneous turbulence driven by an energy input modulated in time is studied 
within a variable range mean-field theory. The response of the system, observed in the second order 
moment of the large-scale velocity difference D(L,t) — ({(u(x + L) ~u(x)) 2 }} oc Re 2 (t), is calculated 
for varying modulation frequencies u) and weak modulation amplitudes. For low frequencies the 
system follows the modulation of the driving with almost constant amplitude, whereas for higher 
driving frequencies the amplitude of the response decreases on average oc 1/uo. In addition, at certain 
frequencies the amplitude of the response either almost vanishes or is strongly enhanced. These 
frequencies are connected with the frequency scale of the energy cascade and multiples thereof. 
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I. INTRODUCTION 

Many turbulent flows are characterized by time depen- 
dent forcing. E.g. the atmosphere of the earth is driven 
by the heating through the radiation from the sun, the 
blood flow in the arteries by the heart beats, etc. Also 
technical flows like the flow in the intake of a combustion 
engine are periodically forced. Another example are es- 
tuaries and adjacent coastal waters, where tidal straining 
leads to a periodic alternation of stratification and tur- 
bulent mixing of saline and fresh water [lj . This results 
in a periodically varying energy dissipation in the upper 
water layers with a 12 hour period. 

The effect of a periodically increasing and decreasing 
energy input on turbulent flow depends on the frequency 
of the driving. This has been studied in reference |2j 
for a turbulent channel flow where the modulations of 
the input rate are generated near the wall. It was found 
that for high frequencies these oscillations are strongly 
damped with distance from the walls, such that they 
do not reach the inner part of the logarithmic boundary 
layer. Another example is Rayleigh-Bcnard convection: 
the interaction between the large scale circulating flow 
and the thermal plumes detaching from the upper and 
the lower boundary layers acts as a stochastically influ- 
enced time-dependent forcing on the turbulent flow in 
the inner region of the cell, as recently shown in |3, 0, HI ■ 
In a von Karman flow between two coaxial corotating 
disks [a, 13 , the energy input rate is not constant if the 
disks are kept rotating at constant speed, but is period- 
ically varying with a geometry-dependent frequency due 
to a coherent vortex precessing around the axis of rota- 
tion. In this case it was also shown, that the statistical 
properties of the turbulent fluctuations are affected by 
the time dependence of the mean flow. However, the av- 
eraged velocity power spectrum still shows Kolmogorov 
scaling over a broad frequency range, in addition to a 
low frequency peak corresponding to the oscillation of 
the mean flow. 

These results raise the question how global quantities 



of a turbulent flow, like e.g. the total energy or the 
Reynolds number, respond to a time dependent energy 
input. This problem is the subject of the present paper. 
From a more fundamental point of view, studying mod- 
ulated turbulence will give more insight into the time 
scales in particular of the turbulent energy cascade. 

In a previous study |8|, the time evolution of the 
Reynolds number in a periodically kicked flow was an- 
alyzed. If the kicking strength and the kicking frequency 
are large enough, the Reynolds number grows and satu- 
rates on a level, which depends on the frequency and the 
kicking strength. The theoretical results from || have 
later been verified numerically in reference [jj- 

In this present paper, we study a related type of forc- 
ing. Rather than periodically kicking the boundary con- 
ditions of homogeneous, isotropic turbulence as in |8|], we 
force the flow through a time-dependent modulation of 
the energy input rate ei n (t) on the outer length scale L, 



i(t) = e (l + esin(wi)). 



(1) 



This means that the flow is stationarily stirred (oc eo) 
to maintain the turbulent flow and, in addition, a time- 
dependent modulation of the forcing (oc e$e) is applied, 
< e < 1. The response of the system to the time- 
dependent stirring can be observed e.g. in the second 
order velocity structure function of the flow field, in par- 
ticular at the outer scale L, D(L,t) = {{{u(x + L,t) — 
u(x, t)) 2 )). This D(L, t) is equivalent to a Reynolds num- 
ber, which W6 d6fill6 &S Rc — Ul.rrns 

L/v. Here, Ui, rms (t) 
is the rms of one component of the velocity, varying with 
time t. Then, disregarding correlations on scale L, 



D(L,t) = 2((u 2 ))=6ul rms (t) = 



&v 2 Re{tf 
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The energy put into the system at time t will travel 
down the energy cascade towards smaller scales and will, 
on average, be dissipated at time t + T, i.e., with a mean 
time delay r. In other words, the dissipation at time t 
depends on how much energy has been in the large scales 
at time t — r. We approximately describe the relevant 
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time scale r for the cascade process by the large eddy 
turnover time tl at that time t — t, 
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,(t-r) y/D(L,t~T)/6' 



(3) 



More accurately, the time scale of the energy cascade 
is given by the sum over the eddy turnover times on 
all decay steps, r ~ X)n T «- I n this sum i the largest 
contribution is the largest eddy turnover time T£. For 
K41 scaling the smaller eddies r n /L — 6 n , where < 
5 < 1, have turnover times r„ = tl5 2tl / z . Thus r = 
tl J2 n d 2n ^ 3 = T L<i- The common choice 6 = 1/2 implies 
a ~ 2.7. Putting into intermittency corrections gives 
slightly smaller values of a. In this present paper we 
shall discuss the influence of a by comparing the limiting 
cases a = 2.7 and a — I. Experimentally, in principle 
the parameter a could be measured by analyzing the po- 
sitions, heights and widths of the response maxima, thus 
giving information about the energy cascade time. 

If the external modulation period a-> _1 is much larger 
than this intrinsic time scale r, lot <C 1, the turbulent 
flow will have time to adjust and will follow the periodic 
variations of the stirring. If, on the other hand, is 
decreased and becomes much smaller than r, the system 
can follow less and less, and feels, at small scales, an 
average time-independent energy input. 

We calculate the time dependence of the response 
D(L, t) — Do(L) to a periodically modulated energy input 
rate, Eq. |JI|I. within a variable scale mean-field theory 
|ldj | for various driving frequencies u>. Here, Dq(L) is the 
second order structure function for a stationary energy 
input rate eo. In general, the energy flow rate through the 
system is an intermittently fluctuating quantity. There- 
fore, the cascade time as well as the response of the sys- 
tem are fluctuating. These fluctuations are neglected by 
the mean-field theory in the present study. However, 
on average these fluctuations result in a mean downscale 
transport of energy which controls the overall properties 
of the flow. Therefore, we believe that within this mean- 
field approach we can grasp the main features of the flow 
correctly. 

The method is explained in the next section. The be- 
havior of the response as a function of the driving fre- 
quency to in the case of weak modulations of the energy 
input rate is analyzed in section II I II In section IIVI we 
discuss an alternative way to introduce time dependence 
into the system. The slightly different case of a modu- 
lated driving force instead of a modulated energy input 
rate is presented in section We summarize our results 
in section EH 



II. METHOD AND MODEL 

In reference [lCj an energy balance equation for the sec- 
ond order velocity structure function D(r) = ({(u(x+r) — 
u(x)) 2 )) for stationary, homogeneous, and isotropic tur- 
bulence has been derived within a variable range mean- 



field theory. Here, u is the velocity and the brackets ((...)) 
denote the ensemble average. One of the essentials of this 
theory is to divide the velocity field into a (spatially aver- 
aged) superscale velocity and a (strongly fluctuating) 
subscale velocity vS r '. The spatial average is performed 
over a sphere of variable radius r, and will be denoted as 
uW(x) = (u(x + y)) { y r) ee ^ f |y| < r d 3 y u(x + y). 

The energy input rate e, n , which in the statistically 
stationary situation equals the total energy dissipation 
rate e, is balanced in accordance with the super- and sub- 
scale decomposition by the energy dissipation rate on all 
scales larger than r complemented by the energy transfer 
across scale r from the super- to the subscales of r. In 
a simplified version the derived energy balance equation 
reads: 



D{rf \ 1 d 



b 3 r 



dr 



D(r), 



(4) 



where v is the kinematic viscosity and b the Kolmogorov 
constant. In the viscous subrange (VSR), where r is 
smaller than the Kolmogorov length scale -q, r < r], the 
dissipation term, i.e., the first term on the rhs of Eq. 
is dominating, and therefore the solution of Eq. J2J is 
D(r) — er 2 /2>v. In the inertial subrange, instead, where 
T) <C r <C L, most of the energy of the eddies is transfered 
down-scale. This energy transfer rate Et, which is given 
by the second term on the rhs of Eq. Q , is determined by 
the decorrelation rate T(r) of the subscale eddies, which 
itself is mainly governed by the energy dissipation rate e, 
see [l0| for details. Note again that in the stationary case 
the energy dissipation rate equals the energy input rate, 
e = e in . In the ISR the second term on the rhs is the lead- 
ing one. Then the solution of Eq. Q is D{r) = 6(er) 2 / 3 . 
The full energy rate balance equation J3J interpolates be- 
tween these two limits. The Kolmogorov constant b can 
be calculated within this theory to be b = 6.3 which is 
consistent with the experimental value [ill H2I IT3L Il4| . 

In our case the flow is not stationary but experiences 
a modulated energy input rate ej n (i). Therefore, e™, 
the structure function D(r), and the dissipation rate e in 
Eq. will depend on time. Furthermore, an additional 
term on the rhs of Eq. (JlJ appears, taking into account 
the non-stationarity of the flow: 
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(5) 



The correlation of the superscale velocities can be writ- 
ten as {{u^{x,t)-u^ r \x,t))) = {{u 2 (x,t))) - + 

2/2, i))t/i'')j/2' ) - Following the arguments in for the 
derivation of Eq. we neglect multiple spatial av- 

eraging, i.e., ({D( yi + y 2 , t)) ~ {D(y,t))i r) . 

In the stationary case the energy dissipation rate e = 
^^f^ff 1 )) can ^ e related to the large scale quantities by 
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D{Lf/ 2 
&I 2 L ' 



(6) 



Extending this expression to the time-dependent case, we 
have to take into account that the energy which is fed into 
the system on large scales at a time t will be dissipated 
on small scales at a later time t + r. We model this as 
follows: The energy dissipation rate at time t is assumed 
to depend on the large scale quantities at time t — r: 



e(t) = c e (D(L,t - t)) 



D(L,t-r) 
6 3 / 2 L 



3/2 



(7) 



c e is a dimensionless function which is approximately con- 
stant (p± 1) for very large Reynolds numbers |15l llq . In 
[l7L Il8l Hi^ l it was shown that in general c e depends on 
the Reynolds number, and therefore on D(L). We here 
use an approximation of the expression derived in |17| for 
high Reynolds numbers: 




The delay time r is determined by the implicit time-delay 
equation ©. Assuming that the solution of Eq. in 
the ISR, D(r) — &(er) 2 / 3 , is valid up to r = L, we can 

write D(r) — (f;) 2 ^ 3 D{L). Within our model, where 
we connect small and large scale quantities at different 
times, the structure function on scale r < L at time t will 
depend on the large scale structure function at an earlier 

time t — T, i.e., we introduce D{r, t) — (f-)"^ 3 D(L, t — r) 
into Eq. (J5J. After multiplying with r, Eq. (JSJ can be 
integrated from r — up to the outer length scale r = L: 

Id,,, , . . D(L,t-r) 3 / 2 

--(D(L,t)-aD(L,t- T) ) = - Lb3/2 



3vD(L,t- t) 
2L 2 



(9) 



where a = || originates from the integration. In 10] it 
has been shown that, in the isotropic and homogenous 
independent of the scale r as the forcing is 
assumed to act on the largest scale L only. In the sta- 
tionary case the lhs of Eq. (O vanishes, and together 
with Eqs. Q and |(SJ|, Eq. JjJJ corresponds to e = e in . 
Eq. 10 contains only large scale quantities. Effects of 
fluctuations in the energy input rate on the statistical 
properties of the turbulent flow as observed in [6| would 
influence the scaling behavior of D(r, t) on intermediate 
scales r and therefore lead to different values of the factor 
a, but the structure of Eq. @ would remain the same. 

Using Eq. J2J), we express the second order structure 
function D(L,t) in Eq. 10 in terms of the Reynolds 



number Re(t): 
L 2 d 



v dt 



(Re 2 (t)-aRe 2 (t-T)) 



2 /6 



3/2 



3 \b 

-6Re 2 (t-T) 



{Re 2 {t-T)f 2 
(10) 



2 e Q L 4 , 
+- — =-(1 + esinwt). 
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Here, we have inserted the time-dependent energy input 
rate, Eq. Q . In the case a of constant energy input rate, 
i.e., e = 0, Eq. simplifies to 







2 /6\ 3/2 o 2 2L 

-(-) Rel-QRel + -- 



(11) 



relating the stationary Reynolds number Reo to the sta- 
tionary input rate, -p-eo(i?eo) = c e (Reo)ReQ. Introduc- 
ing the reduced Reynolds number R(t) = Re(t)/Reo and 
the non-dimensional time i/ r L as ^ (analogously for r 
and ui), Eq. 110(1 becomes 



d(R 2 {t)~aR 2 {t-T)) 
dt 



2 /6 



3 V b 



3/2 



(i? 2 (t-T)) 3 / 2 



Re, 



-R 2 (t-r) 



(12) 



+ (l 6 / & ) 3/2+ i)( 1+esin ^- 

Here, r° = — is the large eddy turnover time of the 

stationary flow. R(t) is of order one. The delay time r 
in units of the time scale r£ is given by 



(13) 



R{t - t) 



Eq. 1(12(1 describes the time evolution of R 2 (t), which is 
the square of the Reynolds number of a flow exposed to 
a modulated energy input rate (Eq. Q), normalized by 
the square of the Reynolds number of a flow where only 
a constant, time-independent, forcing is applied. 



III. RESPONSE OF TURBULENT FLOW TO 
ENERGY INPUT RATE MODULATIONS 

A. General trend 

In the present study we shall restrict ourselves to the 
case of weak amplitude modulation, i.e., e in Eq. is 
small. Then we expect that also the oscillating response 



A(t) = R 2 (t) - 1 



(14) 



has a small amplitude, and we can linearize Eq. 1(12(1 . 
The time delay r is approximated by a time-independent 
constant which in our time units is simply a. This 
approximation is justified as long as |A| -c 1. In section 
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IIII CI we shall discuss the limits of this approximation. 
We first consider a = 1 which means that the cascade 
time r is taken as the large eddy turnover time r£. The 
resulting equation of motion for the response A(i), 



dt 



(A(f)-aA(t-r)) 




r) 



esin tot, 



(15) 



can be solved analytically. The solution to the linear 
equation l|15|l can be calculated using the ansatz: 



A(t) = eA{uj) sin (cot + . 



(16) 



Here, ^4(w) is the amplitude, and </) is the phase shift 
of the response which also depends on to. Inserting this 
expression into Eq. I|15|l gives the explicit solution of the 
linear response equation l|15|) : 



A(i) 



\3\bJ ^ Re 



— cos u>t + a cos 



U7 + fit 



sin 



+ r) 



(1 



or 



(!) 



6\3/ 2 



2a coswt 



,(!) 



(17) 



In the following, we set the Kolmogorov constant 6 = 6 
for simplicity, which is near to the calculated value 6.3 
Q0| and to the experimental value in the range 6 — 9 
To recover the expressions for a gen- 
eral b one has to replace in the following results the 
terms (1 + and (| + ^) by ((6/b)^ + ^) and 

(|(6/6) 3 / 2 + -j^), respectively. The mean amplitude 
of the response is determined by the energy input rate 
(| + -^)e, i.e., the last term on the rhs of Eq. fT5|l . The 
time derivative on the lhs of Eq. I|15|l leads to a mean de- 
crease of the amplitude as 1/lo. Due to the two terms in 
Eq. (|15JI containing the time delay r = a, corresponding 
terms in the second fraction of the solution l|17|) appear, 
cx a and cx (1 + -^-)/cj, respectively, which, by the peri- 
odic dependence on lot induce a periodic variation of the 
amplitude with the frequency to. For low frequencies the 
terms oc (1 + 6/Reo)/u, originating from the first term 
on the rhs of Eq. (|15fl . dominate, whereas for high fre- 
quencies the terms oc a, due to the second term on the 
lhs of Eq. Ijl5|) . become more important. The latter, in 
particular, lead to a periodic variation of the response 
amplitude up to very high frequencies. 

The linear response A(t) oc e of the flow (with Re^ = 
10 4 ) is plotted in FigQ] for four different modulation 
frequencies. Also the modulation of the energy input 
rate, ei n (t)/eo — 1 is plotted in FigQ The deviation 
of the Reynolds number from its stationary value Reg, 
A(t) = (Re 2 (t) - Refy/Rel, oscillates with the same fre- 
quency as the driving, for all frequencies to. The am- 
plitude A of this oscillation depends on the frequency. 
For the two small modulation frequencies, u> = 1CP 3 and 
to = 10 , the amplitude of the response A(t) is nearly 
the same, about two thirds of the amplitude e of the 
driving. For higher frequencies, the amplitude A of the 
response decreases. In the case of to = 10 we observe 
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FIG. 1: Response A(t) (solid lines) for four different mod- 
ulation frequencies ui, the time dependent part of the en- 
ergy input rate, ei n (t)/eo — 1 (dotted lines). The modulation 
amplitude is 10% of the constant input rate, e = 0.1, and 
the Reynolds number of the stationary system is chosen as 
Re = 10 4 . (a) lutI = 10" 3 , (b) lut^ = 0.1, (c) lot° l = 10, (d) 
lot° l = 100. 



a phase shift between the forcing and the resulting re- 
sponse. FigEJi shows the amplitude A{lo) as a function 
of the driving frequency for Reo = 10 4 . For low frequen- 
cies the amplitude remains constant, and is two thirds, 
whereas for large frequencies the amplitude of the re- 
sponse A(i) decreases oc 1/iu. In addition to this decrease 
we note certain frequencies for which the response ampli- 
tude becomes large or very small. The distance between 
two maxima or two minima of the amplitude is nearly 
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FIG. 2: (a) Amplitude A of the response A(t) as a function of 
the driving frequency u (log-log-scale) for weak modulations 
(e = 0.1) of the input rate ei n , and Reo = 10 4 . The time 
scale of the energy cascade is set to r = a = 1. The dashed 
line denotes the low frequency limit of the oscillation ampli- 
tude, 2/3, and the dotted line corresponds to the mean trend 



of the high frequency limit, 4;- 



Inset: linear-scale-plot of 



the response amplitude versus frequency. The small arrows 
indicate the frequencies u) r (in units of t£) of the response ex- 
trema calculated from the extrema of the denominator in Eq. 
I)17fl . The horizontal arrow denotes the frequency distance 
Alo (in units of r£) between two frequencies for which the 
amplitude is maximal (or minimal). It is Aui ~ 2n/r for high 
frequencies, (b) Same as (a) but with a cascade time scale 
r = a = 2.7 different from the large eddy turnover time t£. 
Note the shift of the response maxima, the less pronounced 
height and greater width of the first, and the more pronounced 
second response peak. 



modulation time scale is much larger than the large eddy 
turnover time, 1/lo 3> 1, i.e., if the driving frequency is 
very small, then the solution fTTj) can be approximated 
by 



(18) 



A(t) ~ e-sin ( u(t + r 



We conclude A — 2/3, while the phase (j> = lot is linear 
in lo for small frequencies. 

If, on the other hand, the modulation frequency be- 
comes very large, i.e, the time scale of the driving is much 
smaller than 1, we see from Eq. (|17fl that the amplitude 
of A decreases as oc 1/lo: 



A(i) 



(1 + _6_ 

^ Ren 



COS Lot 



(oo(t + t) 



1 + a 2 — 2a cos lot 



(19) 



^ 3 + Ren > 



■^j of this high frequency 



The mean trend oc 

limit is also plotted in FigEt- The crossover between the 
regimes of Eq. (|18|l and l|19|) takes place at Lo cross ~ 1. 
This can be seen in FigEK- The crossover frequency is 
not changed by taking into account the cascade time r = 
a 7^ 1 , as can be seen in FiglSb which shows the response 
amplitude as a function of frequency for a = 2.7. 

We have considered here only the case, where the Kol- 
mogorov constant 6 = 6. For a general 6, the crossover 
frequency is at Lo cross ~ (6/6) 3 ^ 2 , as can be seen from the 
solution Q17[). This means, that the crossover from the 
regime of constant amplitude to the regime of 1 /w-decay 
takes place at a smaller frequency if b is larger. The posi- 
tions of the response maxima, however, are only slightly 
shifted by a different b. 

In conclusion, as long as the modulation frequency of 
the energy input rate is smaller than 1, i.e., the large 
eddy turnover time is shorter than the period of the forc- 
ing, the system has time to follow the periodic modu- 
lations with an almost constant amplitude. For higher 
frequencies instead, the oscillations become too fast for 
the system to follow, and therefore, the response becomes 
weaker and weaker, and phase shifted. Then the system 
experiences the fast modulation more and more as a con- 
stant average energy input, and the oscillations of the 
response vanish as 1/lo. This high frequency behavior 
has also been found for spin systems driven by an oscil- 
lating magnetic field [20| . 



constant, see the inset of Fig|2k. This periodic behav- 
ior in the w-dependence of the response amplitude is due 
to the time delay r. We shall explain this in the next 
section. 

There are three time scales in the solution (|17|) of Eq. 
()15JI : The large eddy turnover time, by definition 1, the 
time delay r — a, which represents the cascade time, and 
the time scale of the external modulation 1/lo. If the 



B. Response maxima 

In FigEJwe have seen that there are certain frequencies 
for which the amplitude of the response becomes large or 
very small. Mathematically, these response extrema orig- 
inate from the minima and maxima of the denominator 
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inEq. |T7|), 




1 + a 



.» -Ren 
-LOL COS LOT — 1 SlllCJT 



(20) 



We calculate the extrema of N(u) numerically. The first 
few of them are indicated by the small arrows in FigEJi. 
The lowest frequency is near to uj r \ ~ 7r/3r ~ 1. There, 
the first and strongest maximum of the response can 
be observed, where the amplitude becomes as high as 
A ~ 4.2. Note, that this frequency is nearly equal to 
the crossover frequency uj cross between the low and high 
frequency regimes of Eq. 1)18(1 and 1(19(1 only in this par- 
ticular case, where a = 1. If we assume an energy cas- 
cade time t = a = 2.7 the frequencies of the maxima are 
shifted towards smaller frequencies. The height of the 
first maximum is decreased, i.e., A ~ 1.2, whereas the 
height of the following maxima is slightly increased, see 
FigEp- For very large frequencies, u > 1, we can esti- 
mate the frequencies of the response extrema also analyt- 

ically. Then the two terms in the denominator cx — 



can be neglected, and the extrema of N(u>) can be ap- 
proximated by the extrema of coswr, 



u! r (n) 



7T 
71- , 

r 



n = 0,±l,±2, . 



(21) 



Now the amplitude of A is at maximum for frequencies 
oj r (n) with even n, and at minimum for io r (n) with odd 
n. The distance between two maximum (or minimum) 
amplitudes is 2tt/t as indicated by the horizontal arrow 
in the inset of FigEl F° r the first maxima and minima at 
moderate frequencies this estimate is an approximation 
only; also their distances are not yet constant as they are 
for high frequencies. 

In the high frequency limit, the oscillation of the re- 
sponse at the frequencies ia r of maximum or minimum 
amplitude is phase shifted by (j) r {m) = (2m + 1)tt/2, 
m = ±1, ±3, 



A(t) 



(I + 7^") (-1 ±a)cosuj r t 



(1 -F a ) 



oc sin (oj r t + (fir)- 



(22) 

The prefactor (—1 ± a) is always negative, i.e., at the re- 
sponse extrema we have A(t) cx — cosuj r t — sin(uj r t+<fi r ). 
In Fig|3|the phase shift <fi{u>), calculated from the solu- 
tion 1(17(1 . is shown as a function of the driving frequency 
uj for all frequencies. As the phase shift starts with 
<fi(u> = 0) = and changes continuously with increas- 
ing frequency, we find that only m = 1 is possible for 
the phase shift <fi r at the response extrema. The frequen- 
cies of the maximum and minimum amplitudes of A are 
indicated by arrows. The only exception is the first max- 
imum, where the approximation for u> r , Eq. I(21|) does 
not yet hold. There, the phase shift is near to ir/2, cor- 
responding to 7ti — 0. Another phase shift in this model 
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FIG. 3: (a) Phase shift 4>(u>) as a function of the modulation 
frequency ui for weak modulation strength e = 0.1, and Reo = 
10 4 . The upper (lower) arrows indicate the frequencies of 
maximum (minimum) amplitude of the response. For small 
u> the phase <fi(u) oc ujt behaves linearly. 



is the one between the response A(t) and the energy dis- 
sipation rate e(t). According to Eq. (J7J) the dissipation 
rate is phase shifted by — lot with respect to the response 
A(t), i.e., this shift is linearly growing with increasing 
frequency to. At the response maxima and minima the 
phase shift is — lo t t ~ — nir. 

The physics behind these response extrema can be ex- 
plained as follows: The time delay r can be regarded as 
the (average) time which the input energy needs before 
it is dissipated at small scales. In the case of maximum 
amplitude of the response the time delay r is a multiple 
jT of the period T = 2tt/lj of the forcing, whereas for the 
frequencies of minimum amplitude the delay r has an ad- 
ditional T/2. Therefore, at the extrema of the response, 
the energy dissipation rate and the response are either in 
phase (maxima) or anti-phased (minima). In the latter 
case the oscillation of the response is strongly reduced. 
If, on the other hand, the driving frequency is such that 
the response and the dissipation rate are in phase, the 
transport of energy through the system is very effective 
and leads to an enhanced oscillation. At the response 
maxima as well as at the minima the phase shift between 
energy input rate and response is <j) r = 3ir/2. 



C. Quality of the approximation for the delay r 

In the above calculations we made an approximation 
for the time scale r of the cascade process. In the lin- 
earized model, we assumed r to be constant, r = tq = a. 
Now we check a posteriori the quality of this approxima- 
tion. The solution l(17|) of the linearized equation 1(15(1 is 
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used to compute the "correct" delay time r step by step: 
The next approximation for r is 



n(t) 



(23) 



where the delay in Eq. (|13fl is still neglected. Further 
steps are: 



r 2 (t) = 
r 3 (t) = 



Vi + ACt-n)' 

a 

Vl + A(t-r 2 j' 



(24) 



etc. 



In Fig^ t q = a, t%, t 2 , and r 3 are plotted for differ- 
ent frequencies. For uo — 0.01, the difference between 
Ti, T2, and T3 is not visible. The variation of the Ti(i), 
(i = 1,2,3), is largest at the frequency where the am- 
plitude of A is maximum, i.e., at ui ~ For all other 
frequencies, including at the response maxima, the vari- 
ation of the Ti(t) is much smaller than To and 1/uj. At 
these frequencies it seems reasonable to approximate r by 
the constant To = a. In Eq. Ijl5(l the delay r enters into 
two terms, in oc <9 t A(t — r) on the lhs, and in oc A(t — r) 
on the rhs. We calculate the relative error of these terms 
if t = To instead of r = tj (i — 1, 2, 3) is employed, using 
the solution {TTJ) for A: 



J Q 27r ^[cos (w(i — r ) + <j>) —cos (uj(t — Ti) + (j))] 2 dt 



for the term on the lhs, and 



(25) 



\ 



/o 



sin 



(w(t- t ) + </!>) -sin (w(t-Ti) + </))] 2 ^ 



p2Tv/u -2(1, \ 
Jo Sm ^(*-7D)- 



dt 



(26) 

for the term on the rhs. The errors Si and 82 are summa- 
rized in table [I] for the four chosen frequencies of FigQ] 
As expected, the errors are largest for the frequency with 
maximum response amplitude, uj — 1.06, where it be- 
comes up to 23%. For oj — 10 and beyond it is between 
1 and 5%. If one would allow for a time dependence of 
t in Eq. (|15f) the response maxima would probably be- 
come broader, possibly less pronounced. However, within 
this mean field theory we anyhow can make only approx- 
imate statements about the frequencies and the values 
of the amplitudes at the response maxima. Namely, in 
the mean field approach the effects of the fluctuations 
on the structure function are neglected. Therefore we 
believe that even with this approximation for the delay 
time t we can qualitatively predict the basic features of 
the system, which are the decrease of the amplitude of 
the response for high modulation frequencies, and the 
existence of response maxima at certain frequencies due 
to the finite time needed by the energy cascade process. 
The validity of the approximation for r will improve for 
smaller amplitudes e of the modulation. However, for 
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FIG. 4: Successive approximations of the delay time t: First, 
constant approximation to = a (solid lines); second, time de- 
pendent approximation n (dotted lines) ; third approximation 
T2 (dashed lines); fourth approximation T3 (dashed dotted 
lines) for the delay time r, see Eqs. H23I24H . (a) tori = 0.01. 
(b) ujtI = 1.06. (c) ljtI = 10. (d) ujtI = 10 2 . In (a),(c) and 
(d) the time dependent Ti(t) for i = 1,2,3 are indistinguish- 
able. 





Si(n) 


Si(t 2 ) 


Si(r 3 ) 


ujtI = 0.01 


2.9 x 10 -4 


2.9 x 10 -4 


2.9 x 10~ 4 


lot'I = 1.06 


0.15 


0.22 


0.23 


urt = 10 


0.016 


0.016 


0.016 


lutI = 100 


0.046 


0.046 


0.046 




«a(n) 


<$2(t2) 


S 2 (r 3 ) 


ujtI = 0.01 


1.7 x 10 -4 


1.7 x 10~ 4 


1.7 x 10 -4 


ujtI = 1.06 


0.22 


0.11 


0.12 


lutI = 10 


0.013 


0.013 


0.013 


utI = 100 


0.036 


0.036 


0.036 



TABLE I: Relative errors Si, 82 according to Eqs. 125|l . l|26|l 
made in the two relevant terms of Eq. 1151 by using the 
constant time delay to = a instead of the higher order ap- 
proximations Ti(t) for T. 



smaller e the total amplitude eA of the response will de- 
crease as well and finally the amplitude of the response 
maxima and minima will become so small that, in exper- 
iments or numerical simulations, the fluctuations will be 
larger than the maxima and minima. 



IV. AN ALTERNATIVE ARGUMENT TO 
INTRODUCE THE TIME-DELAY 

The energy balance equation J3J and the expression for 
the energy dissipation rate e, Eq. © , hold for stationary 
systems. The time dependence of the quantities in these 
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equations in section [H] has been introduced a posteriori 
by arguments based on the picture of the energy cascade. 
It was not derived from the Navier Stokes equation, but 
is a modeling ansatz. Therefore, there are several argu- 
ments to introduce this time dependence. We want to 
discuss here another way of arguments which leads to a 
slightly different equation for the response. The idea is to 
start from an equation which is already integrated over 
all scales, i.e., does not depend on the scale r any more 
in contrast to Eq. 10} in section [n] The total energy per 
unit mass of the flow is E ~ ?>u\ rms /2. It is basically 
determined by the energy of the large scales. The change 
with time of this energy equals the dissipation rate and 
the energy input rate: 



dt 



E{t)=-e(t)+e in (t). 



(27) 



As the energy needs a time r to travel down the eddy 
cascade before it is dissipated, e at time t may be ex- 
pressed with Eq. and E = \D{L) by the total energy 
E at time t — r: 



5(t)=c e (E(t-r)) 



3/2 



E(t - r ) 
L 



3/2 



Together with the approximation for c £ , Eq. (JSJ, we get: 



d F m- {AE{t 



)) 



3/2 



-6—E(t-r) 



■e in (t). (28) 



As in section [H] we express the energy E by the Reynolds 
number, E — |^-i?e 2 , write the energy input in terms 
of the stationary Reynolds number Re^, Eq. (|llfl . 
and introduce the reduced Reynolds number, R(t) — 
Re(t) / Reo. Then, in time units of t l : 



3/2 



(R 2 (t-r))^-^-R\t-r) 
Re 



esinwt). (29) 




The only difference between this equation and the previ- 
ous one, derived in section [D] (Eq. flUO , ^ s that here the 
term oc dR 2 (t — r)/dt is missing. 

If we solve Eq. (|29|l within the same linear approxima- 
tion as employed in section ITTTI for Eq. (|12|l . we find the 
same features for the response, see Fig. OD The solution 
of the linearized equation obtained from Eq. (|29|l reads: 



A(t) 



(1 + _6\ 



1 + 



coswiH p. sin (w(t + r) 



1 + (3±35l) 9 



■ sin lot 



( 3 °) 

Here, we have again set 6 = 6 for simplicity. The re- 
sponse maxima are also observed, but they are less pro- 
nounced and slightly shifted. The amplitude at the first 
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FIG. 5: Amplitude A of the response A(t) as a function of 
the driving frequency u (log-log-scale) for weak modulations 
(e = 0.1) of the input rate e™, and Reo = 10 4 calculated from 
Eq. 12911 in linear approximation. The dashed line denotes 
the low frequency limit of the oscillation amplitude, 2/3, and 
the dotted line corresponds to the mean trend of the high 
frequency limit, 2 u . Inset: linear-scale-plot of the response 

amplitude versus frequency. The small arrows indicate the 
frequencies of maximum amplitude uj r (in units of t£) calcu- 
lated from the minima of the denominator in Eq. I30H . The 
horizontal arrow denotes the frequency distance Aui (in units 
of t£) between two frequencies for which the amplitude is 
maximal (or minimal). It is Aui ~ 2tv/t for high frequencies. 



(and strongest) maximum has only a value of Ae — 1-6. 
In the linear response solution l|17|) of the previous model 
the terms originating from the second term on the lhs of 
Eq. I|15fl were responsible for the strong variation of the 
amplitude at high frequencies. These terms are missing 
in the present model. Therefore, we observe weaker am- 
plitude maxima and minima at high frequencies in this 
model, cf. Fig|SJ If we take the extended cascade time 
t = a > 1 into acount, e.g. a = 2.7, the response max- 
ima are shifted towards smaller frequencies as discussed 
in section lTlI ti\ However, in this model, the heights of all 
maxima including the first one is then slightly increased. 
At the response maxima the energy cascade time scale 
r and the period of the driving modulation are not mul- 
tiples of each other as they are in the previous model, 
i.e., the response and the energy dissipation rate are not 
exactly in phase. If one would observe the response max- 
ima in experiments or numerical simulations, one could 
distinguish between the two models by studying the ra- 
tio between the frequencies of the response maxima and 
the cascade time scale r. The phase shift <f> between the 
energy input rate and the response becomes negative and 
oscillates around — 7r/2 for higher frequencies. At the re- 
sponse extrema it is near to 4> r ~ —n/2. Note that in 
the previous model the phase shift was always positive. 
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The two arguments to introduce the time delay are 
similar and are based on the same physical idea of a finite 
time lapse of the cascade process. However, we tend to 
prefer the first one, section [H] because it introduces the 
time dependence at an earlier stage. Eq. @ still resolves 
the scales r and it is therefore closer to the Navier Stokes 
equation than Eq. 1|27|) . 

V. RESPONSE OF TURBULENT FLOW TO A 
MODULATED DRIVING FORCE 



Here, /o is the strength of the (stationary) forcing and 
ef the amplitude of the modulation. As has been dis- 
cussed in section [H] we express the response in terms 
of the Reynolds number Re(t) and relate the station- 
ary Reynolds number Reo with the stationary forcing 
strength /o, similar to Eq. Then we introduce 

the reduced Reynolds number R(t) = R ^ and the di- 

mensionless time t = t/r®. The tilde is dropped in the 
following. The analogous equation to l(T2l becomes: 



In the previous sections we have studied the effect of a 
modulated energy input rate on turbulent flow. However, 
the energy input rate may not be a quantity which can be 
easily controlled in experiments. In some experiments it 
is more convenient to modulate the driving force instead. 
Then the resulting energy input rate as well as the total 
energy of the system can be considered as a response of 
the system. Therefore, in this section, we show how to 
treat this slightly modified case within the variable range 
mean-field theory and what differences we expect in these 
two different response functions. 

The derivation of Eq. JUJl for the response of the sys- 
tem in terms of the structure function D{L,t) remains 
the same as explained in section [HI The energy in- 
put rate ei n (t) in that equation is given by ei n (t) — 

{(u'f\x,t)f^ r \x,t))). To introduce a modulated forcing 
instead of a modulated energy input rate, we therefore 
assume: 



,(t) * D(L,t) l / 2 f(t) 

= D(L,t) 1/2 f Q (l + e/sinwi) 



(31) 



d{R 2 (t)-aR 2 {t-T)) 
Jt 



(32) 



^-{R 2 {t-r)fl 2 
(^ 2 W) 1/2 ^ + ^)(l + e/sin^), 



where b is set to b = 6. We again assume small modu- 
lation amplitudes, i.e., e/ <C 1, and linearize Eq. I|32[l in 
A(t) = R 2 (t) — 1. As before, the time delay r is approx- 
imated by the time-independent constant a. With the 
same ansatz Eq. (|16f) as in section UlI Al for modulated 
energy input rate, the linearized equation can be solved 
analytically, and the solution reads: 



A(t) - e, 



(1 + _6_) 

V 3 T Re > 



cos ut + a cos (u(t + T)j + 1 -^ 5SL sin (oj{t + rfj - <3H 2 " e ° ) smut 



1-1 S_ (2 , 

1 + O? + (it^L)> + 



— 2a cos lot — 2 



■ sm LOT 



-) 

— (a sin lot j 22 - cos lot) 

1 (33) 



This solution is very similar to the solution l|17(l for a 
modulated energy input rate, but it contains some ad- 
ditional terms in both the numerator and the denomi- 
nator. These terms only slightly modify the frequency 
dependence of the response A(t). In FigEJthe amplitude 
A[lo) of the response A is plotted as a function of driving 
frequency for Reo — 10 4 . As for the modulated energy 
input rate we note that the amplitude remains constant 
for low frequencies and decreases as oc 1/lo for high fre- 
quencies. Also the response maxima and minima can be 
observed. Quantitatively, the low frequency limit for a 
modulated forcing is different from the modulated energy 
input rate case. For low driving frequencies, lo <C 1/t, 



we can approximate Eq. (|33() by: 
-(1 sin (w(t + r)) 



A(i) ~ ej 



2 

3uj v 



3oj 



sin tot) 



2 + & 



2-^- COS COT 



(34) 



The terms 6 / Reo 1 have been omitted here for simplic- 
ity. In the limit lot —> 0, with sinojr — -> and cos lot — > 1, 
the amplitude A of the response (cf. Eq. I|16|0 becomes 
equal to one instead of two thirds (cf. Eq. 1)181) ) for 
a modulated energy input rate. The frequencies of the 
rcponse maxima and minima are determined by the ex- 
trema of the denominator of the solution (|33|l . They are 
slightly shifted as compared to the case with modulated 
energy input rate (Eq. H17|) ~). The amplitude at the first 
maximum is smaller than in the case with modulated en- 



10 




imated by 



,(*) 



1 ~ ±A(i) + e f smut + 0(A 2 ). The 



-2 -1 

iog 10 (^ L °) 

FIG. 6: Amplitude A of the response A(t) as a function of 
the driving frequency ui (log-log-scale) for weak modulations 
(e/ = 0.1) of the driving force /, and Reo = 10 4 . The time 
scale of the energy cascade is set to t = a = 1. The dashed 
line denotes the low frequency limit of the oscillation ampli- 
tude, 1, and the dotted line corresponds to the mean trend of 
the high frequency limit, f^r- The dashed-dotted line rep- 
resents the amplitude of the resulting energy input rate ei n (t). 
Inset: linear-scale-plot of the response amplitude (solid line) 
and the energy input amplitude (dashed-dotted line) versus 
frequency. The small arrows indicate the frequencies ui r (in 
units of t°) of the response extrema calculated from the ex- 
trema of the denominator in Eq. 1331 . The horizontal arrow 
denotes the frequency distance Auj (in units of r£) between 
two frequencies for which the amplitude is maximal (or min- 
imal). It is Au; ~ 2tt/t for high frequencies. 



ergy input rate, namely Ae — 2.7. However, in the limit 
of very high driving frequencies, Eq. (|33[1 can be 

approximated by Eq. ljT$|) . i.e., the response amplitudes 
of both cases become identical. 

It was pointed out in the beginning of this section that, 
if we modulate the driving force, the energy input rate is 
not a controlled quantity, but can be considered as well 
as a response of the system. This has been measured in a 
recent experimental study by Cadot et al. [2l| . Within the 
mean-field theory the energy input rate for a modulated 
driving force can be calculated as: 



6in,0 



= + A(i)(l + e/ sinwt), 



(35) 



1/2 

where etn.o = ((D^ Q fo)) is the stationary energy input 
rate for constant forcing without modulation. In order 
to extract the amplitude of the energy input rate, we fit it 
by a function of the form e e * = l + efA ein sm(uj(t + (f))). 
This is justified as long as the modulation amplitude 
is small, i.e., e/ <C 1, because then A(t) is of the same 
order of magnitude as ej and Eq. i|35|l can be approx- 



amplitude A ein of the energy input rate is included in 
FigEl as a dashed-dotted line. For low driving frequen- 
cies, u> <C 7) the amplitude A ein is nearly constant and 
is 3/2, whereas for high frequencies it decreases and fi- 
nally saturates at one. Also the response maxima can be 
observed in the energy input rate: At the same frequen- 
cies, where the response A shows amplitude maxima, we 
observe a maximum directly followed by a minimum in 
the amplitude of the energy input rate. 

In conclusion, if the driving force instead of the energy 
input rate is modulated, the general behavior of the re- 
sponse in terms of the second order structure function on 
the larges scale remains the same, including the response 
maxima and minima. For low driving frequencies, the 
amplitude of the response becomes equal to the ampli- 
tude of the forcing. In addition, the energy input rate 
can be regarded as a different measure for the response 
of the system, which also shows the response maxima at 
frequencies connected with the energy cascade time scale 



VI. CONCLUSIONS 

We calculated the response of isotropic and homoge- 
neous turbulence to a weak modulation of the energy in- 
put rate ei n within a mean-field theory. For low frequen- 
cies the system follows the input rate modulation whereas 
for high frequencies the amplitude of the response de- 
creases oc 1/cj. Due to the intrinsic time scale of the 
system, the eddy-turnover time r, which also character- 
izes the energy transport time down the eddy cascade, 
there are certain frequencies, cu r ~ n^, where the am- 
plitude of the response is either increased or decreased. 
At these frequencies the phase shift <fi between the en- 
ergy input rate and the response is <j) r ~ 37r/2. The 
response extrema occur when the eddy-turnover time is 
an even or odd multiple of half the modulation period 
T/2 = tt/lu, respectively. In the case of response max- 
ima, the energy dissipation rate and the response of the 
system are in phase. This can be understood as a very 
effective transport of energy through the system. At the 
amplitude minima, instead, the response of the system is 
strongly reduced. Then, the energy dissipation rate and 
the response are exactly anti-phased. 

In the mean-field approach the fluctuations of the en- 
ergy flow rate through the system and of the large eddy 
turnover time are neglected. In experiment or numerical 
simulation the fluctuations are however present. They 
may lead to broader and less pronounced response max- 
ima, i.e., partly wash out the response maxima and min- 
ima. 

With increasing modulation amplitude e of the energy 
input rate the response maxima are expected to become 
more significant due to the better signal to fluctuation 
ratio. But remember that for higher modulation ampli- 
tudes e, the time scale of the eddy cascade, which enters 
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into our model as a time delay, becomes time dependent. 
This as well could lead to less pronounced response max- 
ima as discussed in section IlII CI 

A way to check if the characteristic feature of the re- 
sponse maxima and minima can still be well identified 
under the influence of fluctuations, would be to perform 
numerical simulations of the Navier Stokes Equation with 
a modulated driving. However, as not only high Reynolds 
numbers are needed to achieve fully developed, isotropic 
and homogeneous turbulence, but also the response as 
a function of time for a wide range of frequencies has 
to be calculated, the computational effort would be too 
high. Therefore, numerical simulations within two dy- 
namical cascade models of turbulence, the GOY shell 

model m EE EE m El and the reduced wave 

vector set approximation (REWA) p^. l30l l3lT ] , were per- 
formed [i^]. These models take into account the fluctu- 
ations. The basic trend of the frequency dependence of 
the response amplitude as calculated within the mean- 
field model can be reproduced in both numerical models. 
We also clearly find the main maximum in both models 
although it os of course washed out by the fluctuations. 
The higher maxima and minima however seem to be com- 
pletely washed out. 

Also a recent experimental study of modulated tur- 
bulence by Cadot et al. [2l| showed evidence for the 
existence of the response maxima. This experiment may 
be comparable with our study of a modulated driving 
force as discussed in section ]V\ The response maxima 



were measured in the amplitude of the energy input rate. 
In addition, a constant response amplitude for low driv- 
ing frequencies and a 1 /w-decay of the velocity response 
amplitude for large frequencies has been observed. This 
is in agreement with the 1/w-decay of the energy re- 
sponse amplitude which we have found in the mean-field 
model. The velocity response (u(t) — uq)/uo — A„(t), 
where u(t) is the measured velocity modulus and uq the 
(stationary) mean velocity, is connected to the energy re- 
sponse A(i) which we have calculated in this paper by 
1+A(t) - u{tf/ul = (1+A u (t)) 2 ~ l+2Au(t)+0(Aj). 
As only small modulation amplitudes are considered the 
term +0(AJ) will be negligible because A u -C 1. There- 
fore, the experimentally measured 1/w-decay of the am- 
plitude A M of the velocity response is just what one would 
expect from our theoretical prediction A oc 1/w for the 
amplitude of the energy response. 

We hope that the present work will stimulate even 
more experimental and numercial studies on the role of 
the energy cascade time scale in modulated turbulence. 
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